Chiral Gap and Collective Excitations in Monolayer Graphene 
from Strong Coupling Expansion of Lattice Gauge Theory 
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Using the strong coupling expansion of the compact and non-compact U(l) lattice gauge theories 
for monolayer graphene, we show analytically that fermion bandgap and pseudo Nambu-Goldstone 
exciton (7r-exciton) are dynamically generated due to chiral symmetry breaking. The mechanism is 
similar to the generation of quark mass and pion excitation in quantum chromodynamics (QCD). 
We derive a formula for the 7r-exciton analogous to the Gell-Mann-Oakes-Renner (GOR) relation 
in QCD. Experimental confirmation of the GOR relation on a suspended monolayer graphene would 
be a clear evidence of chiral symmetry breaking. 
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Graphene is a monoatomic layer of carbon atoms with 
a honeycomb lattice structure^. A novel feature of 
graphene is that electrons and holes at low energy have a 
linear dispersion relation around two independent "Dirac 
points" in the momentum spaced. Then the charge car- 
riers on graphene can be described by massless Dirac 
quasiparticles^. The system has however a critical dif- 
ference from relativistic electrons^; the Fermi velocity v F 
of the electrons on graphene is about 300 times smaller 
than the speed of light c. This leads to an effective 
enhancement of the Coulomb interaction among Dirac 
quasiparticles. In such a strong-coupling situation, elec- 
trons and holes on graphene may form an exciton con- 
densate and create a gap in the fermion spectrum lead- 
ing to semimetal-insulator transition. To show unam- 
biguously that undoped monolayer graphene suspended 
in vacuum becomes an insulator is one of the impor- 
tant theoretical challenges^. Also, this problem has 
much in common with the dynamical breaking of chiral 
symmetry in strongly-coupled relativistic field theories 
such as quantum chromodynamics (QCD) 6 ". So far, var- 
ious theoretical methods such as the Schwinger-Dyson 
equation^—, renormalization group equations^ and lat- 
tice Monte Carlo methodsii^ have been applied to study 
the dynamical formation of the fermion gap in low-energy 
effective theories of graphene. 

The purpose of this Rapid Communication is to shed 
lights on the strong coupling regime of graphene at 
zero temperature from an analytic method of strong 
coupling expansion (See Refs llall4l for its recent ap- 
plications in QCD). We start with a "braneworld" 
or "reduced QED" model of graphene 8,10 in which 
(2+l)-dimcnsional Dirac fermions are coupled to (3+1)- 
dimensional Coulomb field. After discretizing this model 
on a square lattice with compact U(l) gauge field and 
staggered fermio n 11 ' 12 , we carry out an expansion by the 
inverse Coulomb coupling and derive an effective action 
for the fermions. The exciton condensate (chiral con- 
densate) and the fermion gap (chiral gap) are obtained 
analytically from the resultant effective action. Prop- 
erties of the pseudo Nambu-Goldstone (NG) excitation 
associated with the exciton condensation are also stud- 



ied: In particular, we derive a mass formula similar to the 
Gell-Mann-Oakes-Renner relation (GOR) in QCD±£. 

Low-energy Euclidean action for Dirac quasiparticles 
on a graphene interacting with the U(l) gauge field is 
described by 8 -^ 

S E = J2j dx(3) *Pf(D{A A ]+m)iP f 
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where the natural unit (H — c = 1) is taken. Since this is 
an effective theory for quasiparticles in the first Brillouin 
zone of the original honeycomb lattice, it has an intrinsic 
momentum cutoff p A < vr/ahc with ai lc = 1.42 A being 
the honeycomb lattice spacing. 

The three-dimensional and four-dimensional Euclidean 
coordinates are denoted by x^ = (t,xi,X2) and x^ — 
(r, x\, X2, Xs), respectively. The Dirac spinor ipf has 
four components corresponding to 2 (the number of 
sublattices) x 2 (the number of Dirac points). The "flavor 
index" / runs from 1 up to the number of flavors N. In 
this Rapid Communication we specifically focus on N — 
2 case, which corresponds to the monolayer graphene 
with the up and down spins of electrons. The Dirac oper- 
ator is defined as -D[Ai] = 74(^4 + ^4) + v F (-fi 91+72^2), 
where A4 is a temporal component of the gauge field. 
The gauge coupling constant for the suspended graphene 
is g 2 — e 2 /e with e being the electric charge and e be- 
ing the vacuum permittivity, g 2 is reduced by the factor 
2/(1 + e) on a substrate with e being the dielectric con- 
stant of the substrate^. The Hermitian 7 matrices obey 
the standard relation {7^,71/} = 2<5 M „. The Fermi veloc- 
ity reads v F = (3/2)ta\ lc = 3.02 x 10~ 3 in the unit of 
light velocity, with the hopping parameter t ~ 2.8 eVjil 
obtained from the spectral slope observed on substrate. 
The bare mass m corresponds to an explicit bandgap 
which may be formed artificially on epitaxially grown 
graphene on substrate^ or on graphene nanoribbon and 
nanomesb^. 

Due to the small Fermi velocity, electron interactions 
are dominated by the Coulomb interaction, so that the 



2 



spatial components of the gauge field, Aj=x,2,a, can be 
neglected. With scaled variables, r — > t/v f and A 4 — > 
v F A 4 , Dirac particles have an effective mass m* = m/v F 
and an effective coupling g 2 — g 2 /v F which is about 300 
times larger than the Coulomb coupling strength in the 
vacuum. In the chiral limit (m — > 0), Eq.(fT]) is invariant 
under U(4) chiral transformation with 16 generators 2 ^,: 
(1, <r) <g> (1,73,75,7375) with 75 = 74717273- Absence of 
73 in -DLA4] is the reason for such large chiral symmetry. 

A regularized form of Eq.(fT]) on a hypothetical square 
lattice with a lattice spacing a = ir/p A reads 12 
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Here all the dimensionful quantities are scaled by a. The 
U(l) gauge action So is written in terms of a time-like 
link variable U±(x) = exp(i6(x)) with — n < 9 < it. The 
fermionic action Sf is written in terms of the staggered 
fermion X through 

M(x) = J2xb(x) X b(x), (4) 

b 

V^(x) = ^2Tj fl (x)xb(x)U tl {x)xb(x + fi), 

b 

V-{x) = E^( a; )Xfc(^ + A)^(^)Xfc(^), (5) 



with fi = 1,2,4 and Ui,2(x) = 1. 6 is the staggered fla- 
vor index which runs from 1 to N/2, since 2 3 doublers 
emerging from one staggered fermion on a 3-dimensional 
square lattice can be identified with four Dirac compo- 
nents times two "flavors" n i 12 . The monolayer graphene 
corresponds to N = 2. The staggered phase factors t]^ 
are 774(2:) = 1,171(2;) = (-1) T , 172(2:) = {-l) T+x \ and 
173(2;) = (-1)^1+^ = e ( x y 

In the chiral limit, fermion action for each flavor 
in Sf is invariant under U(l) v xU(l) A chiral trans- 
formations; {xb{x),Xb{x)) ^> (e ,? vxb(x),e~^v^ & (x)) 

and (xb(x),Xb(x)) A (e l ^( x ) Xb ( x ) , e l ^( x ) x b (x)) . 
These are remnants of global U(4) chiral symmetry of 
Eq.JT])£L Under the U(1) A rotation, we have M(x) — > 
e Ziti K <*)M{x) and V±(x) -> V±(x), so that the chiral 
condensate (xx) serves as an order parameter for the 
spontaneous symmetry breaking, U(l) v xU(l) A — >U(1) V . 

We define an expansion parameter, /3 = 1/g 2 , so that 
the strong coupling limit corresponds to f3 — > 0. (By tak- 
ing v F ~ 0.003, (3 in the vacuum is estimated as 0.04, 
while that on the SiC>2 substrate is 0.1.) Then, the par- 
tition function expanded by Sq ~ 0(/3) becomes 
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FIG. 1: Induced four-fermion interaction in the strong cou- 
pling expansion. The open (filled) circle represent % (x)- ( a ) 
In the LO, the time-like links (red arrows) in Sf cancel with 
each other to leave a spatially local interaction, (b) In the 
NLO, the time-link links in Sf are canceled by those in So 
(blue arrows) to leave a spatially non-local interaction. 



The 8 integration can be analytically performed order by 
order in /3 for general N—: When the link variables e 
cancel with each other, the fermion self-interactions up 
to 2A r -fermi term are induced. Hereafter, we will focus 
only on the monolayer case (N — 2)^ Up to 0(0), we 
obtain 



s* = E 



\ E (V?(x)-Vf(x))+mJU(a 
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+f E E (V?(x)Vf( x + 4) + (V+ o Vf)) . (7) 

The second line in Eq.© is the leading-order (LO) term 
of 0(/3 ) which is local (non-local) in space (time) as 
shown in FigOJa). The third line is the next-to-leading- 
order (NLO) term of 0(/3) which is non-local in both 
space and time as shown in FigQJb). Note that the 
gauge field propagating along the third spatial dimension 
2:3 starts to appear at 0(/3 3 ) in the strong coupling ex- 
pansion. The non-local four-fermi interactions in Eq.([7]) 
can be liniarized by the extended Stratonovich-Hubbard 
transformation 14 : exp(aAB) ~ J[d<pd<p*] exp[— cv(|</?| 2 — 
Atp — Bip*)], where A and B are fermion bilinears and 
a is a positive constant. By introducing two complex 
auxiliary fields 4>(x) and \(x) corresponding to the LO 
and NLO terms and integrating out the fermion fields, 
we arrive at Z = J[d(f>d(f>*}[dXdX*]e- Se!t ^' X \ where the 
axial U(1) A rotation induces the transformation, 4>(x) — > 

In the mean-field approximation where fluctuations of 
4> and A are neglected, free energy per unit space-time 
lattice cell at zero temperature, F c s(4>), can be obtained 
after eliminating A by using the stationary condition 
(5S eS (<p,\)/6\ = 0): 
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FIG. 2: The free energy F e g (<f>) in the lattice unit as a function 
of \cj>\ for 13 = and m = 0. 



Here G 1 (fc; 0) = X) J= i 2 sm2 % + l m * — ^/2| 2 is the t w0 
dimensional bosonic propagator with an effective mass, 

to* — 4>/2 and J k = ^ j_^,^dk\ fj^.^dk^. Alternative 
way to derive Eq.© is to treat the 0(/3) term of Eq.([7]) 
as a first order perturbation. 

The free energy F c g((j)) in the chiral and srtong cou- 
pling limit (to = 0, /3 = 0) is shown in FigJ5] for illus- 
tration. From Eq.©, we find that F e g((j) 00 ) ~ \4>\ 2 
due to the tree-level term, while F e g((f> — > 0) ~ const. + 
|c6| 2 In |c6| 2 due to the fermion one-loop term. Therefore, 
we can exactly show, in the case when N = 2, that 
dynamical chiral symmetry breaking takes place in the 
strong coupling limit. Since the 0(/3) correction from 
the third term in Eq.® grows as \<f>\ increases, the chiral 
condensate, a = \{xx)\> ^ s a decreasing function of /3. Up 
to the linear terms in /3 and to, we have 



a ~ (0.240 - 0.297,3 + 0.0239 ma)a 



(9) 



where we recover the lattice spacing a. If we em- 



ploy a 



'he 



1.39 keV as a typical cut- 



off scale of our effective theory, we obtain a 



[(0.680- 0.421/3- 



1.39 m y 
cV 



keV] . Note that our ap- 



proach is limited to the strong coupling regime, so that 
it is inappropriate to extract the critical coupling /3 C for 
semimetal-insulator transition from Eq.(|9]). The total 
fermion mass Mp is a sum of the dynamical mass and 
the bare mass in G(k; </>), 



Mp = {aa 2 /2)(v F /a)+m, 



(10) 



which reduces to (0.523—0.623/3) eV +3.05to for a ~ a^. 

Let us now consider collective excitations: a phase fluc- 
tuation of the order parameter ( "7r-exciton" ) analogous 
to the pion in QCD, and an amplitude fluctuation of the 
order parameter ( "cr-exciton" ) analogous to the c-meson 
in QCD. In terms of the auxiliary fields, the former (the 
latter) corresponds to ^> n (x) (cj) a (x)) in the decomposi- 
tion, 4>(x) = a + 4> a (x) + ie(x)4> 7r (x). Propagators of the 
collective modes within the one-loop approximation in 
the strong coupling limit (/3 = 0) reads 



D ^{P' iuj *) 

1 1 + coshu;* 
~ 2 8 



(11) 

H(k,p; a)G(k;a)G(k+p; a), 



where H(k,p;a) — J2j=i 2 sm sm (fcj + Pj) ± (to* + 
a/2) 2 with the + (— ) sign corresponding to the 7r-exciton 
(the cr-exciton). The 0(/3) correction to the above ex- 
pression can also be obtained. The actual dispersion rela- 
tion without the scale transformation of the time variable 
r is derived from the pole of the Euclidean propagator, 

In the chiral limit, 0^(0,0) = is identical to the 
gap equation, dF cS (a) j da = 0, so that the 7t-exciton is 
indeed a NG boson associated with dynamical breaking 
of chiral symmetry. For the leading order in to, the ir- 
exciton mass, M v = w„.(p = 0), reads, 



J V M™=° V a J 



(12) 



As long as < to < 2 meV is satisfied for a ~ 
ahc, M n < Mp holds, so that the 7r-exciton is the 
lightest mode in the system. The relation, M„ oc 
yfm, is similar to the GOR relation for the pion ob- 
tained from current algebra in QCD&i 5 -: Indeed, the 
axial Ward-Takahashi identity for the present system is 
((^J£ xial (x) - 2mP{x))P{y) - 2M(y)5 xy ) = with the 
axial current J^ xial (x) = ^e(x)(V~(x) - V+(x)) and the 
pseudoscalar density P(x) = ie(x)M(x). Saturating this 
identity by the 7r-exciton and using Eq. (fT2"| j£ we obtain, 
in the leading order of to, 

(F;A^) 2 = ma, F; = aa 2 (8v F a) 



-1/2 



(13) 

where a takes the value in the chiral limit (to = 0) and 
the temporal "pion decay constant" F£ is defined by the 
matrix element, (0| Jf xial |7r) = 2F£u) n . The dispersion re- 
lation for the 7r-exciton is also obtained from Eq.(fTT|): At 
low momentum with /3 = m = 0, we have ~ v„ \p\ with 
the pion velocity v n — 4.69t> F = 0.0141. As for the mass 
of the cr-exciton, we obtain, M a ~ (1.30 — 0.47/3) (v p /a) + 
22.6to by solving D7^(0,iM a /v F ) = 0. This is compa- 
rable to the cutoff energy scale of the present lattice, 
E A = v F (n/a). Thus the result is not universal unlike 
the 7r-exciton case, and analogous to the situation for the 
broad a- meson in QCD 6 . Shown in Fig|3]is a qualitative 
summary of the spectra of the fermion and collective ex- 
citations obtained in this study. 

As shown in Refsfl2l. I23I and [24], one may employ a 



non-compact formulation of the gauge action, S G 

. -,2 



(NC) 



2.3 



in order to avoid 



e(x)-6(x + j) 

anomalous phase transition from magnetic monopole 
condensation- 25 . Within the NLO, we find that the 
non-compact results are obtained by the rescaling, e.g. 
cr( NC )(/3) = cr( c )(2/3) + (9(/3 2 ). Our relations a^ c \(3 = 
0) = cr( c )(/3 = 0) and <j( nc )(/3) < cr (c) (/3) in the strong 
coupling region are consistent with the recent lattice 
Monte Carlo simulations 2 ^. 

In this Rapid Communication, we performed an an- 
alytical study of the monolayer graphene in the strong 
coupling regime of U(l) lattice gauge theory. An ef- 
fective action at zero temperature for the Dirac quasi- 
particles is derived up to next-to-leading order of the 
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FIG. 3: Typical mass gaps for the fermion (Mf), the tt- 
exciton (M T ) and the cr-exciton (M a ). The cutoff energy 
of our effective theory reads E A = v F (ir/a) < v f {-k /ate) ~ 
10 eV. 



evidence for the dynamical chiral symmetry breaking in 
monolayer graphene. 

There are numerous problems to be examined in the 
future. Generalization of our approach with the tadpole 
improvemen t 23 ! 24 can be performed. To study the renor- 
malization effect on v F , we need to consider the excitonic 
fluctuations acting on Dirac quasiparticles. These fluctu- 
ations are also important for thermal phase transition of 
graphene from insulator to semimetal. To study universal 
low-energy behavior of the graphene, it would be useful to 
construct a chiral effective theory for light 7r-excitons, a 
non-covariant analogue of the chiral perturbation theory 
in QCD 26 . To be more faithful to the U(4) chiral sym- 
metry at low energies in Eq. ([1} , we should employ lattice 
gauge theory with domain- wall or overlap fermions. Fi- 
nally, to study multilayer graphene where the inter-layer 
electron hopping depends on how the layers are stacked, 
it would be important to develop a lattice gauge theory 
preserving the original honeycomb structure. 



strong coupling expansion. Dynamical breaking of chi- 
ral symmetry and associated formation of a chiral gap 
are found. We showed that the 7r-exciton (similar to the 
pion in QCD) behaves as a pseudo Nambu-Goldstone 
boson in the strong coupling regime. A mass formula for 
the 7r-exciton analogous to the Gell-Mann-Oakes-Renner 
relation in QCD is derived. If the GOR-type mass for- 
mula of the 7r-exciton can be experimentally confirmed, 
e.g. through transport phenomena, it would be a good 
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